Astron. Nachr. / AN , No. ,1 -g]( ) / DOI please setDOI! 



The black hole in NGC 1313 X-2 
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The amount of data available for NGC 1313 X-2 make it a cornerstone for the study of ultraluminous X-ray sources 
(ULXs). We modelled the optical and X-ray data of this ULX with a binary evolution code that takes into account X-ray 
irradiation. We restricted the candidate binary system to be either a 50-IOOA/q black hole (BH) accreting from a 12- 
15Mq main sequence star or a 20-Mq BH with a 12-15Mq giant donor. If the orbital period of the system is ~6 days, a 
20Mq BH model becomes unlikely and we are left with the only possibility that the compact accretor in NGC 13 13 X-2 
is a massive BH of 50— lOOiW©. We briefly discuss these results within the framework of an alternative scenario for the 
formation of ULXs, in which a portion of them may contain BHs of <; 30 — 40Mq formed from very massive stars in a 
low metallicity environment. 
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1 Introduction 

Several pieces of observational evidence strongly suggest 
that a large fraction of ULXs are accreting black hole X-ray 
binaries with massive donors (see, e.g., Zampieri & Roberts 
120091 and references therein). Their X-ray spectral proper- 
ties, along with the long-term flux variability, are similar 
to those observed in Galatic X-ray binaries although, in the 
last few years, Chandra and XMM-Newton observations re- 
vealed new behaviours, showing the existence of spectral 
states with rather peculiar properties (e.g. Gladstone et al. 
2009). At the same time, while most ULXs appear to show 
little variations on timescales of seconds to hours (e.g. Heil 
et al. 120091 1, in some cases fast X-ray variability is rem- 
iniscent of the behaviour observed in Galactic accretion- 
powered sources (rms variability, broad band noise, quasi 
periodic oscillations; e.g. Strohmayer & Mushotzky 2003 



Strohmayer et al. 120071 Feng et al. 120101 Heil & Vaughan 
120101) . Probably, the most compelling evidence of the bi- 
nary nature of ULXs is provided so far by the detection of 
a modulation in the X-ray light curve of M 82 X-l (~ 62 
days; Kaaret et al. l2006l and NGC 5408 X-l (~ 115 days; 
Strohmaver l20 1 Oi l, that are interpreted as the orbital periods 
of these systems. 

A number of ULXs are associated to stellar optical coun- 
terparts and some of them have been identified with stars of 
known spectral type (e.g. Liu et al. 2002 2004, Kaaret et al. 
120041 Mucciarelli et al . 120051 Soria et al. 120051 The coun- 
terparts appear almost ubiquitously hosted in young stellar 
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environments (e.g. Pakull et al. 120061 Ramsey et al. 2006 
Liu et al. 120071 ) and have properties consistent with those 
of young, massive stars. NGC 1313 X-2 is one of the few 
ULXs with a well identified optical counterpart. It is located 
~ 6' south to the nucleus in the barred spiral galaxy NGC 
1313 at a distance of 3.7-4.27 Mpc (Tully [19881 Mendez et 
al. 120021 Rizzi et al. 120071 ). Its observed X-ray luminosity 
varies between a few xlO 39 erg s _1 and ~ 10 40 erg s _1 in 
the 0.3-10 keV band (Feng & Kaaret l2006l[2007l l. The large 
amount of data available make this source a cornerstone for 
the study of ULXs. 

2 A bit of history: X-rays 

NGC 1313 X-2 was discovered by Einstein (Fabbiano & 
Trinchieri 1987 ) and included in the Einstein Extended Medium 
Sensitivity Survey as MS 0317.7-6647. Stocke et al. (1995) 
originally proposed that it could be either a Galactic iso- 
lated neutron star or a binary containing a massive BH in 
NGC 1313. 

It was then repeteadly observed with other satellites, in- 
cluding ROSAT (1991-1998; Stocke et al. [19951 Colbert et 
al. [19951 Miller et al. [19981 Schlegel et al. 120001 ASCA 
(1993-1995; Petre etal. [T994l Makishimaet al. [2000b , Chan- 
dra and XMM-Newton (2000-2006; Zampieri et al. [20041 
Feng & Kaaret 2005 and references therein). The X-ray 
spectrum of NGC 1313 X-2 can be modelled in different 
ways and shows the typical features shared by other bright 
ULXs, such as the presence of a soft, thermal component 
and a rollover of the spectrum above ~ 3 keV (e.g. Glad- 
stone et al. 2009 , Stobbart et al. 2006 and references therein). 
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The amount of XMM-Newton spectra now available allow 
for a tentative identification of different spectral states (see 
Pintore & Zampieri, these Proceedings). 



3 A bit of history: optical 

The optical counterpart of NGX 1313 X-2 was first identi- 
fied on an ESO 3.6 m R band image thanks to an accurate 
Chandra astrometry and after registering the X-ray image 
on the position of SN 1978K (Zampieri et al. 120041 ). ESO 
VLT images of the field showed that the counterpart was 
actually composed of two distinct objects (CI and C2), sep- 
arated by ~ 0.7" (Mucciarelli et al. [2005 ). Further refine- 
ment in the astrometry and accurate modelling of the optical 
emission indicated that object CI was the more likely coun- 
terpart (Liu et al. 120071 Mucciarelli et al. 120071 Patruno & 
Zampieri 2008 ). However, this was established beyond any 
doubt by the detection of the 4686 He II emission line in its 
optical spectrum, a characteristic imprint of X-ray irradia- 
tion (Pakull et al. 120061 Grise et al. l2508l . 

The stellar environment of NGC 1313 X-2 has also pro- 
vided interesting constraints. There are two groups of young 
stars spread out over ~ 200 pc. Isochrone fitting of the 
colour-magnitude diagram of these groups has been attempted 
and provides cluster ages of 20±5 Myrs (Pakull et al. 2006 
Ramsey et al. 120061 Liu et al. 120071 Grise et al. I2008T) . As 
several other ULXs, NGC 1313 X-2 is also associated with a 
very extended (~ 400 pc) optical emission nebula that gives 
important information on the energetics and lifetime of the 
system (Pakull & Mirioni 2002). Assuming that it is formed 
in an explosive event or that its mechanical energy comes 
from the ULX wind/jet activity, the characteristic energy 
and age of the nebula turn out to be ~ 10 52 — 10 53 erg and 
~1 Myr, respectively (Pakull et al. 120061) . Finally, several 
different estimates of the mass and luminosity class of the 
optical counterpart of NGC 1313 X-2 were reported in the 
literature (see e.g. Patruno & Zampieri 2010 and references 
therein). 

Recently, Liu et al. (2009) found a possible periodicity 
of 6.12±0.16 days in the B band light curve of the opti- 
cal counterpart of NGC 1313 X-2, that was interpreted as 
the orbital period of the binary system. Three cycles were 
detected in the B band, while no modulation was found in 
V . Previous studies carried out on the available HST and 
VLT observations led to negative results (Grise et al. 2008). 
More recently, lack of significant photometric variability on 
a new sequence of VLT observations has been reported by 
Griset al. (2009). In principle, the detection of the orbital 
period would definitely confirm the identification of the op- 
tical counterpart and the binary nature of this system. Most 
importantly, it would open the way to perform a dynamical 
measurement of the BH mass. 



4 Optical constraints 

In the following we use the X-ray luminosity and all the 
available optical data on NGC 1313 X-2 to constrain the 
properties of this ULX and its BH. This is probably one 
of the most accurate analyses that is possible to perform at 
present on a single ULX. We adopt the V and B band pho- 
tometry of the counterpart as determined by Mucciarelli et 
al. 120071 As the source is variable, we have further corrected 
the magnitudes and colours by using the average value of V 
and by propagating the errors on V and B. The error in the 
absolute magnitudes is taken to be equal to the maximum 
uncertainty in the different distance determinations of NGC 
1313 (Tully [l988l Mendez et al. 120021 Rizzi et al. 120071) . 

Concerning the reddening, two different estimates of the 
colour excess were adopted in the literature: E(B — V) = 
0.1 d2007l Grise et al. 120081) and E{B -V) = 0.3 (Liu et 
al. [2007b . In the following we chose E(B - V) = 0.11. 
If we assume E(B — V) = 0.3, it is not possible to ob- 
tain agreement between the binary evolution models and all 
the simultaneous constraints coming from observations. For 
E(B — V) = 0.11 the adopted absolute V magnitude is 
therefore in the range -4.38 to -4.79 and the B — V colour 
is -0.13 ±0.06. 

Three different values of the metallicity inferred from 
an abundance analysis of the HII regions in NGC 1313, all 
sub-solar, are reported in the literature: Z = 0.1 Zq (Pi- 
lyugin|200B, Z = O.2Z (Ryder [19931 and Z = 0.4Z© 
(Walsh & Roy [19971 Hadfield & Crowther [20071). There- 
fore, although the metallicity in the environment of NGC 
1313 X-2 appears to be definitely sub-solar, the actual value 
is still uncertain (see Ripamonti et al., these Proceedings; 
see also Pintore & Zampieri, these Proceedings, for an in- 
dependent metallicity estimate from X-ray spectra). 

Finally, we re-analyzed the optical photometry of NGC 
1313 X-2 with the aim of clarifying the statistical signifi- 
cance of the orbital periodicity identified by Liu et al. ( 2009 ). 
As reported in Impiombato et al. (these Proceedings), the 
folded B band light curve shows a 6 days periodicity with 
a significance slightly larger than 3o\ This suggests that the 
periodicity may be there, but the low statistical significance 
of the B band modulation, along with the lack of detection 
in the V band, make its identification uncertain. 

5 Modeling the NGC 1313 X-2 binary system 

We modelled NGC 1313 X-2 using the binary evolution 
code described in Patruno & Zampieri (2008), that com- 
putes the evolution of the orbital parameters and the com- 
panion star taking into account the contribution of the ac- 
cretion disc and X-ray irradiation during the contact phases. 
For this analysis, we developed a new version of the code, 
with a more accurate calculation of the optical colours (Pa- 
truno & Zampieri l2010l) . 

The evolution is computed using two different values 
of the (zero age) metallicity, both sub-solar: Z = O.2Z 
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Fig. 1 CM diagram for binaries with a 2OAf BH under- 
going case AB mass transfer. The black triangle marks the 
position of the optical counterpart of NGC 1313 X-2. All the 
tracks are plotted only during the contact phases. The lower 
part of each track (up to the turn around) refers to the MS 
RLOF phase, whereas the upper part to the H-shell burning 
RLOF phase. The labels on each curve indicate that mass of 
the donor star (in M ). 
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Fig. 2 Same as FigureCQfor a 7OAf BH (case AB). 

and Z = 0.5Zq. We found that the effects on the evolu- 
tionary tracks and the stellar colours are minimal. Hence, 
for reference, in the following we will consider only tracks 
computed for Z = O.5Z . 

If the donor starts a contact phase via Roche-lobe over- 
flow (RLOF), we assume that a standard accretion disc forms 
around the BH. The accretion rate is instantaneously taken 
to be equal to the mass-transfer rate from the companion. 
When the accretion rate M exceeds the Eddington rate AlEdd, 
we impose M = M.Edd and assume that the excess mass is 
expelled from the system. In our simulations, all the binaries 
that start RLOF on the main sequence (MS) have also a sec- 
ond episode of mass transfer after the terminal age main se- 
quence (TAMS). Following Patruno & Zampieri (2010jl, we 
term such binary models as Case AB, while those starting 
the first contact phase after the TAMS are termed Case B. 



For the latter, typically the contact phase is too short when 
compared to the characteristic age of the nebula (~ 1 Myrs), 
which provides an estimate of the active phase of the ULX. 
Therefore, from now on, we will focus on CaseAB systems. 

5.1 Evolution of a binary with a 20M Q BH 

FigureQ]shows the evolutionary tracks on the color-magnitude 
(CM) diagram computed for binary systems with a BH of 
20Mq and donors of different masses. Considering the in- 
trinsic variability of the counterpart, its position on the CM 
diagram may be broadly in agreement with 10 — 2OM 
donors during H-shell burning. However, only stars in the 
12 — 15Af have characteristic ages consistent with the es- 
timated parent cluster age (15-25 Myrs). Donors this mas- 
sive, during H-shell burning, provide a mass transfer rate 
largely above Eddington. In this case there is excess mass 
that must be somewhat expelled from the system. Beam- 
ing and/or super-Eddington emission is necessary to explain 
the observed ULX luminosity. We note also that the system 
certainly underwent a first contact phase during MS, during 
which it had time to energize the surrounding nebula. 

5.2 Evolution of a binary with a 7OAf BH 

Figure [2] shows the evolutionary tracks on the CM diagram 
for a systems with a BH of 7OAf and donors of differ- 
ent masses. The position of the optical counterpart is now 
consistent with the tracks of stars on the MS, in the mass 
range between 12 and 20Mq. Again, the further constrain 
inferred from the parent cluster age (15-25 Myrs) restricts 
the allowed donor masses to the interval 12 — 15Af . The 
mass transfer of a 15A/ at this stage is typically between 
2 x 1CP 7 and 2 x 10~ 6 A/ yr _1 , consistent with the accre- 
tion rate (~ 1CP 6 A/ yr _1 ) inferred from the measured 
average luminosity of NGC 1313 X-2 (~ 4 x 10 39 erg 
s ; Mucciarelli et al. 2007). Similar results are obtained 
for BHs with a mass between 50 and 100 M . 

5.3 Evolution of the orbital period 

The evolution of the orbital period P is rather different for 
the two systems. In order to have a first contact phase dur- 
ing MS, the initial orbital separation is ~ 1 day. For the 
2OA/ BH binary, the orbital period at the end of MS is — 4 
days. During the H-shell burning phase it increases rapidly 
up to 20 days. When the 12 — 15M donor tracks cross the 
position of the counterpart on the CM diagram, P > 8 days. 

For the 7OM BH binary, the orbital period at TAMS 
reaches ~ 6 days and is between 5.5 and 6.5 days when the 
evolutionary tracks of the 12 — 15M donors are consistent 
with the photometric properties of the optical counterpart. 

6 Discussion 

According to the actual value of the orbital period of the sys- 
tem and considering the constraints from the inferred mass 
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transfer rate, the position on the CM diagram, the character- 
istic ages of the parent stellar cluster and the bubble nebula, 
the following scenarios are possible: 

- P ^ 4 days: mass transfer from a ~ 15M© MS donor 
onto a ~ 2OM BH, if all the mass is accreted onto the 
BH and super-Eddington emission is allowed. 

- P = 5 — 7 days: mass transfer of a MS donor of 12 — 
15M Q onto a 50 - 1OOM BH. 

- P *i 8 — 10 days: mass transfer of a H-shell burning 
donor of 12 - 15M onto a - 2OM BH. If P ~ 12 
days, the photometry is consistent with the possition of 
the counterpart only for donors of ~ 2OM , too young 
with respect to the age of the parent cluster. 

- P 15 — 20 days: mass transfer of a H-shell burn- 
ing donor of ~ 15M onto a ~ 2OM BH, only if X- 
ray irradiation is switched off (because, e.g., emission is 
beamed and not able to hit the disc and donor surfaces). 

If we use the tentative identification of the ~ 6 days 
orbital period originally reported by Liu et al. (2009) as a 
further constraint, the observational data appear to be con- 
sistent only with the second scenario and we are left with 
the only possibility that the compact accretor in NGC 1313 
X-2 is a massive black hole of 50 — 100Af Q . 

Such a system would fit well within a massive BH sce- 
nario for ULXs forming in a low-metallicity natal environ- 
ment (Mapelli et al. 120091 Zampieri & Roberts l2009l Mapelli 
et al. 120 101 Mapelli et al., these Proceedings). Albeit uncer- 
tain, in fact, the metallicity of the environment of NGC 1313 
X-2 is definitely sub-solar. As mass loss becomes less effi- 
cient in such environments, already at Z Ss 0.2 — O.3Z 
massive stars may retain a significant fraction of their en- 
velopes at the end of their life (see e.g. Zampieri & Roberts 
2009 and references therein). If the envelope is more mas- 
sive than ~ 30 — 4OM , a low metallicity star may col- 
lapse directly to form a BH (Heger et al. 1200 3 , Belczynski 
et al. 120091 1. This occurs because, after core collapse, the 
supernova shock wave loses too much energy in trying to 
unbind the envelope until it stalls and most of the star col- 
lapses into a BH of mass comparable to the star final mass 

30 - 4OAf ; Fryer [19991 Zampieri [2002>. With these 
massive BHs only very modest beaming or slight violations 
of the Eddington limit (a factor of a few) are needed to ac- 
count for the luminosity of bright (Jb 10 40 erg s" 1 ) ULXs. 

As noted above, however, the identification of the or- 
bital period of NGC 1313 X-2 remains uncertain. Further- 
more, we did not perform a complete survey of the param- 
eter space evolving case AB systems with BH masses be- 
tween 20 and 5OM . For BH masses in this range there 
may also be agreement with observations with other values 
of P. Finally, modelling the amplitude and shape of the light 
curve (considering X-ray irradiation and ellipsoidal modu- 
lations) may give additional constraints. A systematic in- 
vestigation of this type is postponed to when a more robust 
assessment of the orbital period will be available. 
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